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Abstract 

A Hamiltonian cycle of a graph is a closed path which visits each of the 
vertices once and only once. In this article, Hamiltonian cycles on planar 
random lattices are considered. The generating function for the number of 
Hamiltonian cycles is obtained and its singularity is studied. Relation to 
two-dimensional quantum gravity is discussed. 
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Hamiltonian cycles have often been used to model collapsed polymer 
globules 0. A Hamiltonian cycle of a graph is a closed path which visits 
each of the vertices once and only once. The number of Hamiltonian cycles 
on a graph corresponds to the entropy of polymer system on it in collapsed 
but disordered phase. There can be no polynomial time algorithms to deter- 
mine whether the number is zero or not which work for arbitrary graphs[|2[. 
Even for regular graphs (lattices), the number of Hamiltonian cycles is not 
obtained exactly except for a few well-behaved cases P, [|, [5|, || 

In the present work, I consider the problem of counting the number of 
Hamiltonian cycles on planar random lattices, or planar random fat graphs. 
I obtain the exact generating function for the number. 

Let S n be the set of all planar trivalent fat graphs with n vertices possibly 
with multiple edges and self-loops. Graphs that are isomorphic are identified. 
The set S n is the labeled version of S n , namely, vertices of G G S n are 
labeled as 1, . . . , n and Gi, G2 G S n are considered identical only if a graph 
isomorphism preserves labels. The symmetric group of degree n naturally 
acts on S n by label permutation. The stabilizer subgroup of G is called the 
automorphism group Aut G. 

A Hamiltonian cycle of a labeled graph G G S n is a directed closed 
path (consecutive distinct edges connected at vertices) which visits each of n 
vertices exactly once. Hamiltonian cycles are understood as furnished with 
a direction and a base point. The number of Hamiltonian cycles of G is 
denoted by 7~t(G) because it is independent of the way of labeling. See figure 
|] for an example. 

The quantity I study in this work is 

and the generating function 

00 

F{p) = Y,P nF - ( 2 ) 

n=0 

It is illuminating and is useful to rewrite F n as the number of isomorphism 
classes of the pair (graph, Hamiltonian cycle): 

F n = #({(G, C n )\G G S n , C n : Hamiltonian cycle on G}/ ~), (3) 
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where C") ~ (G 2 , C£ ) if and only if G\ and G 2 are isomorphic (forgetting 
the labels) and the isomorphism maps onto with the direction and 
the initial point preserved. Eq. ([3]) implies that F n is an integer though the 
definition ([J) involves a fraction. 

The equality (|[) can be shown as follows. 

C n )\Ge S n , C n : Hamiltonian cycle on G}/ ~ 

=— : Hamiltonian cycle on G} 

' Ges n 



TJ.I Z ' 



n! 

n! G ^„ #Aut G 



ft! 



W(G). (4) 



The first equality follows from the fact that there are n! inequivalent labeling 
of (G, C). In the last line, the definition of Aut G is made use of. 

To compute F(p), the method used in refs. [0, |[ H is followed. One 
walks along the Hamiltonian cycle in the specified direction starting from the 
base point and records the order of right and left turns. Then one obtains 
a diagram consisting of T's as depicted in the center of figure 0. A cycle 
corresponds to exactly 2 out of 2 n diagrams consisting of n T's. The two are 
mirror images of each other. There are m diagrams that have k openings on 
the left hand side and n — k on the right hand side. To reproduce the graph 
completely, one has to connect the n openings pairwisely. There should be 
no connection between the right hand side and the left hand side because the 
cycle divides the sphere into two disks. The connection pattern should be 
able to be drawn on a disk faithfully, i.e. without intersection. The number 
of ways of contracting k objects on a disk is denoted by A^. Then one has 



n 1 , 

x — > 1 in 



^2\kj AkAn - k - (5) 

k=0 v ' 

The factor 1/2 comes form the identification of the mirror images. There is 
no double count besides it. The number has been obtained by Brezin et. 



al. using the gaussian matrix model [10] : 



1 f 2 , . fO (k : odd), 

A k = — dxWA-x 2 x k = < fe! / 1 . \ ( 6 ) 



2W-2 ' " 1 (A:: even) 
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Figure 1: An example of planar 
trivalent fat graph G with n = 4. A 
Hamiltonian cycle is drawn in thick 
line. For this graph, 7i(G) = 4! 
because I distinguish the directions 
and the base points of the cycle. 



Figure 2: Two of contributing di- 
agrams for n = 6. The sequence 
of T's in the center is obtained by 
walking along the Hamiltonian cycle 
(the thick line in the T's). Openings 
in both sides should be connected by 
arcs. On the left hand side, there is 
a single possibility because A 2 = 1. 
On the right hand side, = 2 and 
the two possible patterns are drawn. 



Plugging this into eq.(p|), one has 



n=0 ^ J-2 71 J -2 

I I f 2 1 f 2 1 

= / dxy A — x 2 — / dy\/A — y 2 — -, 

22tt J_ 2 2W-2 l-A{x,y- P y 



where 



By using 



A(x, y\ p) = px + py 



(7) 



(8) 



1 

1 — px — py 



da e ( - 1+p{ * +2/))c \ 



(9) 
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one obtains 



i r°° -a 
F(p) = — 2 da —{h(2pa))\ (10) 
Z P Jo a 

where I\ is the modified Bessel function of the first order. 

From the representation (^), it is apparent that F(p) diverges at p = 1/4. 
The singularity of F(p) brings information of large- n behavior of F n . By 
studying it one can inspect the properties of large graphs. One can show 
that F(p) has a singularity 

F{p) = ^(4p — l) 2 log \Ap — 11 + regular terms. (11) 

lQnp 6 

by examining (^). This result is in accord with the scaling behavior of a 
polymer loop with a base point on planar random surfaces in the dense 
phase 

F(p)-(p-\? (12) 

obtained in J7J. 

Each planar trivalent fat graph corresponds to a triangulation of a sphere 
under the dual transformation. Thus F(p) is the partition function of two- 
dimensional quantum gravity with an unusual weight. Triangulations that 
do not admit a Hamiltonian cycle are excluded from the path integral. On 
the other hand, ones admitting Hamiltonian cycles are weighted by 7i(G). 
Eq. ([TTD shows that the string susceptibility exponent 7 st ring for this system 
is zero. Because there is a base point on the Hamiltonian cycle, one has a 
local degree of freedom which gives rise to a factor proportional to the area 
of the surface. In this viewpoint, one should compare the value 7 str ing = 
with scaling exponent of c = —2 quantum gravity with a puncture operator 
inserted. 

This system can be generalized by introducing two parameters p\ and P2 

as 

oo ^ 

JWa) = £ £ #AulG^' m(G) ' (13) 

k,m=0 GeS k+m 

where Tik,m{G) denotes the number of Hamiltonian cycle with k right turns 
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and m left turns. Then straightforward calculation shows 

F{pi,p 2 ) =7;7T / dx\/A-x 2 — / dy^A-y 2 - 

2 27T J_ 2 27T J _ 2 l-pix 



PlV 



00 e 



da^ r I 1 (2p 1 a)h(2p 2 a) 



(2pi + 2p 2 — l) 2 log \2pi + 2p 2 — 1| + regular terms. 



-1 c 3/2 3/2 

lo7rp 1 p 2 



(14) 



The function F(pi,p 2 ) is not only symmetric but with a singularity described 
in terms of p\ +p 2 only. This implies that the sum is dominated by the cycles 
with an equal number of right and left turns. The fluctuating geometry does 
not change the dominance of the maximum of the factor ( fc ^ m ) at k — m. 
I comment that in flTB| ) the integer ^ir(k — m) can be considered as the 
holonomy || of the cycle if one regards the system as a quantum gravity. 

This solution can be extended to planar g-valent fat graph. One should 
just replace A(x,y;p) in (|7|) by A(x, y;p) = p YH=o x s y q ~ 2 ~ s - One can even 
consider the system with 

JV 

M?, y\ {Pi,A) = Pij^y 3 - ( 15 ) 

i,j=0 

The system now describes Hamiltonian cycles on graphs consisting of (i + 
j + 2)-valent vertices with various The weight pij corresponds to the 
2 71 " x .^^ 2 -turn on a (i + j + 2)-valent vertex. For example, if one takes 
A(x,y;pi,p2) = p\x 2 +P2xy + piy 2 , Hamiltonian cycles with bending energy 
parameterized by P1/P2 are realized on a random 4-valent lattice. 

When one generalizes the system as ([15|), the critical line of F({pi t j}) is 

A{x,y,{ Pi>j }) = l. (16) 

Assume that there exists q — 2 e N such that 

p itj = 0tti + j>q Q -2, (17) 
p ijj ^0ifi + j = q -2. (18) 

This amounts to assuming that the largest coordination number of the graph 
is go and the cycle can turn arbitrarily on the vertex. From direct evaluation 



5 



2 




A(x,y;(p))=l 



Figure 3: The shaded region is the domain of integration (|7|) . The singularity 
occurs when the line px+py = 1 touches this shaded region. For the general- 
ized case (|I5D , the line is replaced with an algebraic curve 1 = A(x,y, {pi,j}) 
but the local situation around (x, y) = (2, 2) is the same as the original case. 



of singular part of ([F|) with (|15|), one can show that the singularity of F({pi j}) 
is still of the form 

F({ Pij })~(5p) 2 \og5p, (19) 

where Sp is a function of Pijs which vanishes linearly on the critical curve. 

This fact can be understood as follows. Let us look at the original case (0) 
with (H). The range of integration of x, y is the square domain —2 < x, y < 2. 
The singularity occurs when the line 1 — px — py = touches this domain 
as depicted in Fig. It occurs on (x,y) = (2,2) assuming p > 0. In the 
case (pf), the critical curve is now an algebraic curve (|16D . It intersects the 
square at (x, y) = (2, 2) transversally. The curve can be approximated by a 
straight line with a slope when one is interested in critical behaviors. The 
slope is the only relevant parameter at the criticality. 

In conclusion, I have obtained the generating function for the number of 
Hamiltonian cycles on planar random lattices and have considered the limit 
of large graphs. 

I thank Mitsuhiro Kato for useful discussions. This work is supported by 
CREST from Japanese Science and Technology Corporation. 

While this paper is being typed, a preprint by B. Eynard, E. Guitter, 



and C. Kristjansen has appeared [11]. It contains interesting results closely 
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related to the present paper. Their definition of Hamiltonian cycle is slightly 
different from that taken here; they do not associate directions and base 
points to cycles. 
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